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CONSTRICTIONS 


By CARL E. KINDSVATER,' M. ASCE AND 
ROLLAND W. CARTER,” A. M. ASCE 
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SyYNopsIs 


Area constrictions in many forms are encountered in both natural and 
artificial open channels. The hydraulics of transitions, control sections, and 
bridge piers and piles has been extensively investigated in the field and in the 
laboratory. It is remarkable, therefore, that most of the investigations of 
obstructions in bridge waterways have ignored the fact that the waterway 
itself is a definitive example of a width constriction in the river channel. The 
bridge waterway is studied in the basic research described in this paper. 

The effect of an area constriction on the hydraulics of steady, open-channel 
flow depends on the geometry of the flow boundaries, the rate of flow, and the 
physical properties of the fluid. The resulting flow pattern is not subject to 
a general analytical solution. A practicable solution of the discharge equation 
has been achieved by the application of a systematic experimental investiga- 
tion to an approximate analysis. Empirical data, which extend the solution 
to a wide variety of boundary conditions, are presented in dimensionless, 
graphical form. A significant feature is a detailed description of the external 
and internal characteristics of the flow pattern. 


NoTaTION 


The letter symbols adopted for use in this paper are defined where they 
first appear, in the illustrations or in the text, and are arranged alphabeti- 
cally for convenience of reference in the Appendix. 
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INTRODUCTION 


Purpose of the Investigation —For several years the Geological Survey, 
United States Department of the Interior (USGS), has conducted field and 
office investigations of various ‘indirect’? methods for determining flood dis- 
charge. Among these methods was the ‘“contracted-opening’”’ method. 
Certain unexplained inconsistencies, revealed in these investigations, resulted 
in the authorization of a program of fundamental research on the contracted- 
opening method. The purposes of the research, begun in February 1951, 
were (1) to acquire a critical knowledge of the mechanics of the flow of water 
through width constrictions in open channels, and (2) to develop, if possible, a 
practicable method (based on field surveys of flood marks and channel char- 
acteristics) for computing the discharge through bridge waterways. 

It was immediately recognized that the attainment of the second objective 
would depend entirely on a preliminary success in accomplishing the first ob- 
jective. It was also realized that the flow pattern involved was not subject 
to a general analysis and that the attainment of the first objective would 
depend largely on a carefully controlled experimental investigation of certain 
basic boundary forms. It was anticipated that the laboratory study would 
not only establish the relative influence of the various flow parameters and 
primary geometric variables, but that a general method of evaluating the dis- 
charge function could be developed and subsequently verified by a systematic 
examination of the influence of additional boundary conditions. Thus, 
during the first phase of the investigation, the ultimate objective was sub- 
ordinated to an intensive investigation of the arbitrarily defined fundamental 
problem of tranquil flow through open-channel constrictions. This part of 
the research is believed to have general significance as a basic tool in applied 
fluid mechanics. 

Scope of the Investigation—The flow of water through bridge waterways 
was the definitive problem in the determination of the scope of the research. 
Normally, this example of rapidly varied motion involves steady, turbulent 
flow, entirely in the tranquil range. The submerged parts of the highway em- 
bankments, bridge abutments, and piers comprise an abrupt width constriction 
in the natural channel. 

The variety of natural boundary conditions, considering the adjacent 
channel reaches and the constrictive elements, is almost unlimited. It was 
therefore necessary to restrict the basic investigation to certain elementary 
boundary forms which could be easily described and conveniently reproduced 
in the laboratory. Thus, the constrictions used in the experimental studies 
were limited to the forms made possible by the use of vertical walls of various 
lengths, widths, and corner roundings. Limitations on the relative size and 
location of the constrictions, the range of stream velocities, depths of flow, 
channel configurations, and channel roughnesses were governed by an arbitrary 
definition of the normal, single-opening bridge waterway. 

The effect of piers and piles, multiple openings, and sloped embankments 
were not investigated during this phase of the research. By restricting the 
study to tranquil flows through constrictions simulating bridge waterways, 
certain problems which are associated with shaped transitions in artificial 
channels were also excluded. These exclusions were justified, in part, by the 
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fact that other investigations*:***’ had provided satisfactory solutions to 


9, ae oe Design of Flume and Siphon Transitions,"’ by J. Hinds, Transactions, ASCE, Vol. 
Pp 
4“The Flow of Water in Fiumes,”’ by F. C. Scobey, Technical Bulletin 393, U.S.D.A., Washington, 
D. C., December, 1933. 
5 “Experiments on the Flow of Water Through Contractions in an Open Channel,” by E. W. Lane, 
Transactions, ASCE, Vol. 83, 1919-1920, p. 1149. 
*“High-Velocity Flow in Open Channels: A Symposium,” ibid., Vol. 116, 1951, p. 265. 


7“Channel Transitions and Controls,” by Arthur T. Ippen. Evgineering edited by H. 
Rouse, John Wiley and Sons, Inc., New York, N. Y., 1950, Chapter VIII, p. 496 


many aspects of the problems related to transitions and nittedtve control 
sections. It is of interest to note that the variety of boundary conditions 
considered is equal to, if not greater than, that covered by available informa- 
tion on the orifice, which is the closed-conduit counterpart of the problem in- 
vestigated. 

Review of the Literature.—Information (published prior to 1953) on the 
theory and methods of computing flow through open-channel contrictions is 
both rare and inconclusive. Ivan E. Houk,* M. ASCE, reporting methods of 


“Calculation of Flow in by Ivan Houk, Technical Reports, Miami Conservancy 
District, Dayton, Ohio, Part IV, 1 


computing flood diese’ tied by the Miami (Ohio) Conservancy District, 
was probably the first to describe the contracted-opening method. Mr. 
Houk credited S..M. Woodward for the procedure used by the District. Mr. 
Houk showed that the discharge equation underlying Mr. Woodward’s analysis 
was identical with that ‘which described the flow through a venturi tube— 
that is, a form of the combined energy and continuity equations. The dis- 
charge equation offered by Mr. Houk included a coefficient which was primarily 
a contraction coefficient. 

E. W. Lane,® M. ASCE, is the first known investigator in the United States 
to have studied simple width constrictions in the laboratory. However, his 
work was concerned primarily with higher Froude numbers than are usually 
involved in natural flow under bridges, and the variety of boundary forms 
studied by Mr. Lane was extremely limited. 

D. L. Yarnell, in his research on flow past bridge piers, made an extensive 
study of the literature, which is annotated in his report. Mr. Yarnell’s ex- 


**Bridge Piers as Channel by David L. Yarnell, Technical Bulletin 442, U.8S.D.A., 
Washington, D. C., November, 1934 


periments ignored the effect of the bridge abutments. For his analysis, Mr. 
Yarnell selected 


“the bridge-pier formulas most commonly used in the United States * * * 
d’Aubuisson’s, Nagler’s, Weisbach’s, and Rehbock’s.’”® 


Mr. Yarnell dismissed as unsound the Weisbach formula, which assumes that 
the total flow through a constricted section is the sum of an orifice discharge 
and a weir discharge. This conclusion, which Mr. Houk had previously 
stated, is generally accepted. The d’Aubuisson formula, in its basic form, is 
simply the discharge equation which results from a simultaneous solution of 
the continuity and energy equations. Therefore, the resulting equation is 
essentially the same as that suggested by Mr. Woodward, and should not be 
termed an empirical “formula.” The Nagler and Rehbock equations are 
essentially forms of the discharge equation which have been adapted specifically 
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to the bridge-pier problem. These equations are more logically termed 
empirical formulas. 

In 1944, the USGS began using a procedure for estimating the flow through 
area-constrictions which was based on a form of the d’Aubuisson equation. 
The USGS (as a result of Mr. Houk’s work and of a limited number of field 
observations) used values of the coefficient of discharge that were related to 
the degree of rounding of the upstream edges of the constriction. 

Laboratory Arrangement.—Laboratory tests made as a part of the investi- 
gation were conducted in the hydraulics laboratory of the Georgia Institute 
of Technology, at Atlanta, through arrangements with the Georgia Tech Re- 
search Institute. The principal item of equipment was a steel flume 18 in. 
high and 10 ft wide, with a usable, uniform-flow length of approximately 21 ft. 
Water was supplied from the laboratory’s recirculating system. Discharges 
were measured by use of either a venturi meter or a bend meter, both loeated 
in the supply piping. Water levels were measured by manometer-connected 
floor piezometers, portable static tubes, and point gages. Velocities were 
determined by a small (Pygmy) current meter and pitot tubes. 

Constrictive elements and false walls were constructed of concrete, alumi- 
num, or exterior-grade plywood. Artificial channel roughness was created by 
various arrangements of chain-link fencing fastened to the smooth, concrete 
floor of the flume. Experimentally determined values of Manning’s n corre~- 
sponding to the three degrees of roughness investigated were 0.012 (for the 
smooth floor), 0.027, and 0.036. Irregular channel cross sections were pro- 
duced by placing various arrangements of } in.-thick, cement-asbestos board 
on the floor of the flume. Manning’s n for the board (smooth side) was deter- 
mined to be 0.011 within the range of conditions tested. 


DESCRIPTION OF THE FLow PATTERN 


General.—There is particular reason for including a detailed examination of 
the external and internal flow characteristics in an investigation of the open- 
channel constriction. 

Since the experimental investigation must be limited to certain basic 
boundary conditions, the successful application of the results to practical 
problems will depend, to a large extent, on a critical knowledge of the internal 
mechanism of the flow. Equally important, however, is an understanding of 
the external characteristics of the flow pattern, so that the data required for 
proper application of the results can be determined. There is ample evidence 
that the failure to recognize the critical variables has accounted for many of 
the inconsistencies attributed to the contracted-opening method of computing 
discharge. 

The External Picture—Fig. 1 illustrates certain basic boundary conditions 
and. hydraulic characteristics of flow through simple width constrictions in 
rectangular open channels. The beginning of acceleration of the central 
body of water approaching an open-channel constriction is indicated by a 
slight decrease in water surface level at a distance upstream from the face of 
the constriction approximately equal to the width of the opening. Decelera- 
tion occurs along the outer boundaries, and a separation zone is created in the 
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corner adjacent to the constriction (zone la in Fig. 1). The size of the separa- 
tion zone and the amount by which its surface level differs from the full stag- 
nation level is a function of the upstream channel characteristics and the 
geometry of the constriction. 

Water from the sides enters the opening as a sharply curved and contract- 
ing stream, indicating that it is being accelerated in directions both normal 
and parallel to the streamlines. Just as the longitudinal acceleration and the 
consequent negative pressure gradient is reflected as a drop in the longitudinal 
profile, the normal acceleration accounts for a considerable difference in level 
between the outside filaments of flow and the central*filaments of flow. It 
generally happens that the greater part of the total drop in average surface 
level occurs very near the upper face of the constriction. 

As the water passes through the constriction, the contracted stream reaches 
a minimum width at section 2 which corresponds to the vena contracta in 
orifice flow. The flow is bounded here by an eddying body of water in zone 3 
which marks a second separation zone. Expansion of the live stream begins 
at the vena contracta and ends at the section of normal regimen in the full- 
width channel section 4. The length of channel required for this process 
depends on the boundary configuration, the roughness of the channel, and the 
magnitude of the pertinent flow parameters. Hunter Rouse,’"° M. ASCE, 


10 “Diffusion of Submerged Jets,”” by M. L. Albertson, Y. B. Dai, R. A. Jensen, and Hunter Rouse, 
Transactions, ASCE, Vol. 115, 1950, p. 639. 


working with a two-dimensional orifice discharging into the atmosphere, found 
that the nominal boundary of the expanding jet would have a slope of 1 on 4 
with respect to the longitudinal axis. H. R. Henry’: working with the two- 
dimensional discharge below a sluice gate (wherein the expansion is limited by 
a fixed lower boundary and a free-surface upper boundary) reported a slope 
of approximately 1 on 6. 

The average level of the stream at the vena contracta (section 2) is ap- 
proximately equal to the level in the downstream separation zone (3). The 
level of the fully expanded jet at section 4 is usually higher than that in zone 3. 
The latter relationship is a function of the shape and friction characteristics of 
the channel between these sections. 

For constrictions of considerable length (measured parallel to the flow), 
eddying water fills the space (3b) between the face of the constriction and the 
contracted stream. Under these circumstances, the level of the water in this 
space is frequently a good measure of the average level of the surface at the 
vena contracta. 

It is apparent, therefore, that certain factors which have a dominant in- 
fluence on the flow pattern deserve closer attention. It is of particular im- 
portance to examine the effect of the fixed boundaries, the separation bound- 
aries, and the upper, free-surface boundary on the distribution of velocities 
and surface elevations over the flow pattern. Under similar circumstances, a 
flow net is frequently utilized. In this case, however, the same restrictions 
which prohibit the accomplishment of a general analytic solution also prohibit 
the construction of a flow net. By temporarily ignoring the influence of the 
free surface, however, an approximate flow net can be drawn on the basis of 
an analogy with the classic, two-dimensional orifice. 
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The Orifice Analogy.—Within the specified limits of the investigation, the 
accelerated motion in the vicinity of an open-channel constriction can be con- 
sidered to be primarily two-dimensional in the horizontal plane. Subject to 
obvious theoretical and practical limitations, it follows that the flow pattern 
can be approximated by a flow net drawn for a two-dimensional constriction. 

Fig. 2(b) shows a thin-plate orifice, or slot, located symmetrically in— 
and normal to—parallel boundaries. The fluid flowing is assumed to be in- 
compressible, and it is assumed to fill the space downstream from the orifice. 
The flow net was constructed by approximate, graphical methods. In its 
preparation, it was assumed that velocities in the approach section were uni- 
formly distributed between the boundaries, but that separation would occur in 
the corner upstream from the constriction. These are contradictory assump- 
tions, because the first implies irrotationality (corequisite to the assumption 
of a non-viscous fluid), and the second implies rotationality. Furthermore, 
the separation boundary shown in Fig. 2(6) is an arbitrarily assumed line, the 
exact location of which is indeterminate. These contradictions, which were 
required for the purposes of the analogy, are not unusual. The resulting flow 
net is not expected to be seriously in error for the case of smooth boundaries 
and high values of the Reynolds number. 

It should be emphasized that the purpose of introducing the orifice analogy 

is to facilitate a qualitative—rather than a quantitative—enalysis of the 
phenomenon under investigation. In fact, the main purpose is to make Fig. 
2(b) available for the subsequent analysis of the internal flow picture. Fair 
agreement between the three-dimensional, open-channel flow pattern and the 
analogous two-dimensional pattern are expected only when the former is re- 
stricted to tranquil, turbulent flows in smooth, regular channels, with low to 
moderate degrees of channel contraction. The investigation of conditions 
outside this range are benefited indirectly by the analogy. 
_ The Internal Picture —Two specific occurrences of separation are indicated 
in Figs. 1 and 2(b). In the first instance, separation occurs near the corner 
formed by the conduit boundary and the upstream face of the constriction. 
The separation at this point is consistent with the admission that the flow of 
real fluids is accompanied by the existence of a boundary layer. Thus, the 
normal deceleration tendencies in the corner immediately adjacent to the 
boundary (where the limiting boundary-layer velocity is already zero) result 
in the reverse flow which initiates separation. The separation zone, desig- 
nated as zone Ia in Fig. 2(b), is a region of eddying fluid under an essentially 
constant piezometric head equivalent to the piezometric head along the 
bounding streamline. Thus, the separation prevents the occurrence of a full 
stagnation head in the corner. The separation also influences the piezometric 
head and the velocity distribution over the adjacent flow pattern. The velocity 
discontinuity at the separation boundary creates eddies which result in a loss 
of energy. 

The second instance of separation, shown in Figs. 1 and 2(b), is marked by 
the boundary between the contracted stream and the eddying fluid in the 
corners downstream from the constriction. Unlike the upstream separation 
zone, this zone is fixed by the boundary conditions and the assumption of 
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irrotationality of the live stream. Thus, the jet outline is assumed to be 
identical with the free-surface boundary which would exist if the stream were 
not enclosed by the fluid flowing. Here, also, the piezometric head in the 
separation zone is essentially constant and equal to the piezometric head on 
the bounding streamline. From the properties of the flow net, furthermore, 
it is seen that the piezometric head in the separation zone must be equal to 
the piezometric head across the vena contracta. This conclusion is sub- 
stantiated by the observation that from Newton’s second law, a pressure 
gradient cannot exist across the boundary of the jet where that boundary is a 
straight line, for in that instance there is zero acceleration normal to the flow. 
It also follows from this principle that the pressure in the zone adjacent to the 
downstream side of the constriction (zone 3a in Fig. 2(b)) is essentially equal 
to that in zone 3 opposite the vena contracta, because of the negligible ac- 
celerations which exist between these points. 

Expansion of the Contracted Stream.—A description of the internal flow 
picture is not complete without an explanation of the expansion or diffusion 
and the resulting deceleration of the contracted stream downstream from the 
constriction. A knowledge of this process is of little significance in relation 
to the discharge equation, but is prerequisite to an understanding of the energy 
transformations which occur downstream from the constriction. 

In Fig. 2(b) the idealized boundary. which divides the downstream separa~ 
tion zone from the live stream is not indicative of the disturbed flow which 
ordinarily characterizes this region. The boundary actually exists as a fairly 
smooth, steady surface only when a heavy fluid, such as water, discharges into 
a large space occupied by a relatively light fluid, such as air. When (as is 
assumed in this case) the fluid discharges into an enclosed space occupied by 
a fluid of similar density, the extreme velocity gradient at the separation 
boundary’ results in shear stresses sufficient to initiate violent eddy action. 
The end effect is a condition of instability wherein the kinetic energy of the 
live stream is steadily converted into turbulent energy, which is eventually 
dissipated through viscous shear. 

As shear along the separation boundary results in deceleration of the live 
stream, the width of the jet increases in the direction of flow in accordance 
with the continuity relationship. Eddies created along the boundary result 
in lateral mixing of the two bodies of fluid until eventually the diffused jet 
occupies the full area of the cross section. Just as the fluid in the jet is de- 
celerated, the fluid originally in the separation zone below the constriction is 
accelerated, entrained, and carried downstream by the mixing process. For 
steady flow to be maintained within the fixed boundaries, therefore, equal 
quantities of fluid must be exchanged across the boundary of the expanding 
flow. The circulatory flow which occurs in the separation region is thus a 
consequence of the mixing process and the slight pressure gradient which 
exists along the conduit boundary between the constriction and the expanded 
jet. 

Distribution of Piezometric Head.—The assumption of irrotationality within 
the limiting streamlines shown in Fig. 2(b) permits the application of the 
constant-energy principle to all points in the flow. From the properties of 
the slow net, therefore, an equation can be written which relates the variation 
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in piezometric head to the relative spacing of the streamlines. The piezometric 
head at any point on the surface of an open-channel flow (in which pressure is 
hydrostatically distributed) is equivalent to the elevation of that point. This 
assumption is valid if the accelerations in the vertical direction are relatively 
small. In this case, the only region in which vertical accelerations might 
account for minor departures from hydrostatic pressure distribution is in the 
region of maximum total acceleration, near the upstream corner of the con- 
striction. Thus, the orifice flow net, drawn for purposes of analogy, yields an 
approximate picture of the variation of surface level in the vicinity of a similar 
open-channel constriction. 

Since it is implied (in the assumption of irrotationality) that the fluid is 
frictionless,” the energy loss between any upstream section of a particular 
streamtube and a section at the vena contracta can be ignored. Thus, Ber- 
noulli’s equation takes the form, 


v? 
29 


in which V is the average velocity in the streamtube cross section; h is the 
piezometric head, defined as the sum of the pressure head and the elevation; 
and g is the gravitational acceleration. From the principle of continuity, as 
applied to the two-dimensional streamtube of width An, 


q = V An = a constant 


in which q is the unit discharge. From the properties of the flow net, the 
normal distance between the streamlines, An, is equal to the normal distance 
between the potential lines, As, at the same point in the flow. 

Therefore, 


in which the term to the left of the equal sign is a dimensionless piezometric- 
head ratio and the term on the right is determined at every point in the flow 
pattern by the relative spacing of the streamlines. As the geometry of the 
flow net is fixed by the geometry of the limiting streamlines, the piezometric- 
head ration is independent of the absolute magnitudes of the reference head, 
ho, the reference velocity, V+, and the size of the orifice. The distribution of 
the piezometric head along the center line of the conduit and the upstream 
boundary of the orifice is shown in Fig. 2(a). These curves can be assumed 
to be fair approximations of the variation in the relative surface elevations at 
corresponding locations at the open-channel constriction. The ordinate 
scale in Fig. 2(a) corresponds to Eq. 4. 
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_ An: _ As: (3) 
Combining Eqs. 1 and 3 results in 
As, \? 
29 


ConstRICTIONS IN RECTANGULAR OPEN CHANNELS 


Selection of the Independent Variables.—The flow pattern in the vicinity of 
an open-channel constriction is influenced by the physical properties of the 
fluid (assumed to be water throughout this paper), the rate of flow, and the 
shape and roughness of the constriction and adjacent channel reaches. In 
this case the influence of surface tension and compressibility can be ignored. 
In the region of comparatively rapid acceleration near the constriction, the 
influence of viscosity can be assumed to be relatively insignificant. Thus, 
the influence of fluid weight is presumed to be dominant among the physical 
forces which govern the motion pattern. The critical flow parameter, therefore, 
is the Froude number F expressed as 


in which y is a significant depth. 

Because only tranquil flows have been considered, the Froude number is 
' always less than unity; and because the range of depths is fixed by arbitrarily 
defined, normal, natural-channel conditions, the velocity and also the unit 
discharge are limited to comparatively narrow bounds. 

The most distinctive influences, therefore, are those which constitute the 
boundary conditions. Among these influences are many natural boundary 
characteristics which are not independently expressible in terms of simple 
mathematical quantities. Others, such as roughness and channel shape, 
might be incorporated in a uniform-flow equation, but they are conveniently 
omitted from an analysis of the discharge characteristics of the constriction. 
For purposes of an approximate analysis, therefore, certain boundary condi- 
tions which are peculiar to irregular natural channels are excluded. The 
subsequent examination considers only rectangular, level channels. 

It is required to select the variables which play major, independent roles in 
influencing the motion pattern previously described. As nothing in the 
mechanical process of dimensional analysis will indicate the omission of 
critical independent variables, every possible variable must be examined. 
Certain significant flow characteristics, such as water levels and velocities, 
are omitted because they are dependent on variables which are more con- 
veniently incorporated in a practical solution. Certain other variables, 
which are preferred for describing the analogous orifice pattern shown in Fig. 


2(b), are necessarily replaced by variables better suited to the description of 


an open-channel flow. 

Dimensional Analysis—The variables selected and illustrated in Fig. 1 
are believed to be sufficient to describe the average flow characteristics. The 
total discharge is denoted by Q; B is the width of the channel; 6 denates the 
width of the opening; Ah = h, — hz is the difference in piezometric heads (or 
surface levels, assuming hydrostatic pressure distribution) between the ap- 
proach section (1) and the section of maximum contraction (2). The average 
depth of flow at section 2 is represented by y2 (at least one depth must be in- 
cluded, and y:2 is assumed to be best defined as an average); L is the length of 
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the constriction at the opening, parallel to the flow; and r denotes the radius 
of the rounding of the upstream corners of the constriction. The density p is 
not included as an independent variable because it is contained in the terms 


Ah and g. 
With the selected variables assembled in a functional relationship, 


fi (Q, B, b, Ah, Y2, L, f, 9g) =0 


they are combined by methods of dimensional analysis into a minimum number 
of significant ratios, 


*\by: Vg Ah’ byzVgy:’ B’ b’ b’b/ 


The first ratio in Eq. 66 is proportional to a coefficient of discharge, C, and 
the second is equivalent to F for the contracted stream. The other ratios 
are representative of the principal geometric properti¢s of the boundary. 

If Eq. 66 is written in the more significant form, 


b L fr 


it is seen that, for any given form of boundary, the configuration of the free 
surface is a unique function of the Froude number of the contracted stream. 
That is, C is a constant for any given boundary and value of the Froude 


number. 

The Discharge Equation.—The form of the discharge function indicated by 
Eq. 6c is approximated by combining the one-dimensional energy equation 
and the continuity equation between sections 1 and 2 in Fig. 1. It is appro- 
priate, however, to recognize the two major energy losses which occur between 
the nominal approach section (1) and the section of maximum contraction (2). 
The least of these losses, designated as h,., is the eddy loss caused by the turbu- 
lence engendered in the upstream separation zone. The greater loss is h,, the 
boundary friction loss, which is primarily a function of the channel character- 
istics. 

An additional consequence of recognizing the influence of the channel 
configuration and roughness is the introduction of factors that account for 
the nonuniform velocity distribution in the flow cross section. Thus, letting 
a represent the ratio of the average kinetic energy to the kinetic energy corre- 
sponding to the average velocity, the energy equation takes the form, 

2 2 
thi + he t+ he + hy 
Introducing C., the coefficient of contraction, in the equation of continuity 


results in 
Q By Vi = C.b y2 


| 
| 
| 
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Solving Eq. 7 and 8 simultaneously for the discharge, 


V4; 


C.b 
29 (hi — hy — he hp 


Qe 


2 
The terms h, and a, Le are dependent primarily on the resistance character- 


istics of the approach channel. These terms are relatively independent of the 
variables which describe the accelerated motion through the constriction. In 
any practicable computation procedure, therefore, the factors are best evaluated 
separately for substitution in an equation such as Eq. 9. The terms h, and a2, 
however, are much less dependent on channel shape and roughness than on 
the geometry of the constriction. Thus, these factors are functions of the 
variables contained in Eq. 6. Therefore it follows that Eq. 9 can be written as 


2 
If it is assumed that 
(11) 


in which K, is an energy-loss coeffieient. From Eqs. 6c and 10, wherein C is a 
coefficient of discharge in the form indicated by Eq. 6c, 


b Ye A r) 


In the foregoing description of the flow pattern, it was observed that the 
elevation of the water surface at the downstream side of the constriction 
(zone 3a in Fig. 1) is essentially equivalent to the elevation of the water sur- 
face at the section of maximum contraction (2). Under these conditions, it 
is possible to compute y, as the difference between the average level at 3a and 
the average bottom elevation at 2. In terms of natural channel conditions, 
therefore, y: is approximately the difference between the average level of the 
water surface at the downstream face of the constriction and the average 
bottom elevation at a corresponding section across the opening. This is not 
only a very practicable procedure, but for all conditions (except, perhaps, 
extreme values of L/b) the computed depth, ys, is a close approximation for y2. 
Eq. 10 thus takes the form, 


72 
Q=Cbys 29 (ab ~ hy (13) 


in which Ah = (h; — A;), the difference in water-surface elevations between 
the nominal approach section (1) and the downstream side of the constriction 
(3a). The discharge coefficient remains a function of the parameters contained 
in Eq. 12. 
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Experimental Verification of the Analysis—In the preceding development 
the form of the discharge equation was determined by an approximate energy 
analysis. The significant parameters that govern the coefficient of discharge 
were determined from dimensional analysis of the assumed variables. The 
primary purpose of preliminary laboratory tests which were made in connection 
with this investigation was to substantiate the form of the discharge function 
and to establish the relative importance of each of the parameters. The 
functional relationship for C, expressed by Eq. 12, was assumed to be correct 
only to the extent that it was a dimensionally correct combination of the 
selected variables. It was assumed that failure to include all of the pertinent 
variables would be indicated by the inability to correlate all the experimental 
data by use of the derived parameters. 

The preliminary exploratory tests were made in a laboratory reproduction 
of the basic boundary conditions shown in Fig. 1. These tests not only sub- 
stantiated the selection of variables but demonstrated that the test procedure 
coula be simplified by ignoring the influence of certain variables for all or 
part of the practical range. Thus, of the ratios contained in Eq. 12, it was 
indicated that the width ratio, b/B, is of major importance, while F, L/b, and 
r/b are of lesser but appreciable importance in defining the C-function. The 
ratio y2/b, a shape factor, appeared to have a negligible influence for the 
entire range of conditions which were considered to be practical. 

With the ratio y:/b eliminated, the functional relationship for C in Eq. 12 
contains four ratios. Not more than two of these ratios can be shown con- 
veniently, together with C, in any single graphical illustration of experimental 
results. Because of the necessity of developing a practicable computational 
procedure, the experimental data from the initial tests were examined in detail. 
After much study, a procedure was adopted which involved an approximation 
method for determining the complex interrelationships of the principal vari- 
ables. The method was then substantiated by demonstrating that it would 
yield results of sufficient accuracy for the entire practical range of all the 
variables. 

An Approximate Solution for Simple Constrictions.—The procedure adopted 
involved a base curve of the relationship between C, 6/B, and L/b which was 
first determined for certain arbitrary constant values of F and r/b. To com- 
plement the base curve, average ‘‘adjustment’’ curves were developed from 
systematic investigations of the separate influences of the latter ratios. In 
every case, the adequacy of the average curves was checked throughout a 
reasonable range of all the independent ratios. With these curves derived 
from the experiments, a straightforward procedure was achieved for computing 
the essential characteristics of flow through simple constrictions in level, 
rectangular channels 


CONSTRICTIONS IN IRREGULAR CHANNELS 


The General Problem.—An analysis of discharge characteristics and a 
satisfactory computational procedure have been presented for a definitive, 
simplified example of an open-channel constriction. For this basic case, 
the coefficient of discharge was shown to be a function of five independent 


467-12 


| 
| 
| 
| | 
| 


ratios, one of which was eliminated as a result of a laboratory investigation. 
One of the remaining ratios is a flow parameter—characteristic of all gravity- 
controlled phenomena—the Froude number. The other ratios are descriptive 
of the geometry of the constriction. The geometry of the channel in which 
the constriction is located is involved, to only a minor extent, in the evaluation 
of the friction-loss term in the discharge equation. 

However, it is necessary to extend the analysis to the more complex bound- 
ary conditions involved in the definitive practical problem—the bridge water- 
way. Outstanding among. these complexities are the configuration and rough- 
ness of the channel. It is assumed that the friction loss can be computed by 
use of the Manning formula or a similar open-channel-flow formula. But the 
influence of a rough, irregular channel on the efflux characteristics of the con- 
striction is a matter which cannot be evaluated by analytic means, nor can it 
be approximated by analogy with the common orifice. This, then, is the 
principal objective of an extended analysis. 

Additional Variables in the Extended Analysis.—As specified under the 
heading “Introduction. Scope of the Investigation,” this study is restricted 
to certain basic boundary conditions. The simplest of these, as illustrated in 
Fig. 1, were involved in the successful accomplishment of an approximate 
solution. The extended analysis involves irregular channels of varied rough- 
ness patterns, and three additional conditions relative to the shape and arrange- 
ment of the constriction in the channel. 

The channel in the vicinity of natural-channel constrictions is irregular in 
the sense that it is comprised of a main channel and one or two overbank, flood 
channels. The cross-sectional shapes tested in the laboratory are shown in 
Fig. 3 and the corresponding dimensions are listed in Table 1. As indicated 
in Fig. 3 and Table 1, the relative width occupied by the overbank channel 
was assumed to vary, and the constricted waterway was assumed to include 
various proportions of main channel and flood channel. It was also assumed 


that each of the component sections might possess different roughness char- 


acteristics. Thus, the roughness patterns tested in the laboratory were com- 
posed of various combinations of artificial bottom roughnesses, varying, in 
terms of Manning’s n, from 0.011 to 0.036. It was assumed that the channel 
would have the same shape and roughness characteristics downstream from 
the constrictions as upstream. Exploratory tests indicated that major varia- 
tions in the shape of the downstream channel have no significant effect on the 
discharge characteristics of the constriction. 

Anticipating that the channel shape and roughness characteristics do not 
have to be included in the array of independent variables involved specifically 
in the discharge function, only three new boundary parameters need be added 
to those contained in Eq. 12. These parameters are described as eccentricity, 
angularity, and entrance guide walls. 

Whereas it was previously assumed that the constrictions would be sym- 
metrically located with respect to the longitudinal axis of the channel, in the 
extended analysis it will be assumed that they might be eccentrically located. 
The degree of eccentricity, represented by the ratio e, is defined as the ratio of 
the length of the shorter obstruction to the length of the longer obstruction. 
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In the intitial analysis, it was assumed that the parallel faces of the con- 
strictive element were perpendicular to the parallel boundaries of the channel. 


TABLE. 1.—VALvEs OF THE VARIABLES SHOWN IN Fic. 3. 


y, in feet 


2.50 to 10.00 
5.56 to 9.05 


Natural channel boundaries are not necessarily parallel, and the thread of the 
stream is not necessarily perpendicular to the parallel plane of the constriction. 
Thus, in the extended analysis, an angle of skew or angularity, ¢, is defined as 
the acute angle between the plane of the constriction ana a line perpendicular 
to the thread of the stream. 

The third added variable involves the shape of the upstream corners of 
the constriction. Exploratory tests indicated that rounding the downstream 
corners has an insignificant effect on the coefficient of discharge. In the initial 
analysis the upstream corners were taken to be either sharp or rounded to 
circular arcs. In this analysis it was assumed that the corners might have 
diagonal, plane faces, simulating vertical guide walls. A single variety of 
guide walls (set at 45° as indicated in Fig. 7(b)), to be shown subsequently was 
investigated. Verification of the approximate solution as applied to this de- 
sign should be adequate evidence of the validity of the solution. The ratio 
W/b is a measure of the relative length of the guide wall. 

The additional boundary conditions impose certain qualifications on the 
parameters included in Eq. 12. Thus, areas cannot be computed as simple 
products of depths and widths, velocities cannot be assumed to be uniformly 
distributed across the channel, and depths cannot be referred to a horizontal 
bottom datum, The most disturbing consequence, however, is the changed 
significance of the width ratio, b/B. 

The Channel-Contraction Ratio.—The width ratio, 6/B, was previously 
used as a convenient measure of boundary geometry. For rectangular cross 
sections this ratio is proportional to an area ratio. With a uniform velocity 
distribution across the section, the ratio is also equivalent to a discharge ratio. 
If the ratio ¢/Q is substituted for b/B, however, ¢/Q is seen to possess a greater 
significance. As illustrated in Fig. 1, q is that part of the total discharge Q 
which occupies an area of width 6 in the total cross section upstream from the 
constriction. To extend this definition to the general case, q is defined as the 
normal discharge capacity of a channel having the area characteristics of the 
opening and the roughness characteristics of the approach channel. As 
ordinarily defined, the “contraction ratio” is a measure of the relative con- 
striction imposed on agivenchannel. Therefore, if m is defined as the channel- 
contraction ratio, then 


Figure B, in feet b, in feet e= = ro 
3(a) pre 1.00 to 8.00 0.00 to 1.00 
3(b) 1.34 to 5.23 0.79 to 1.00 0.17 to 0.33 
: 3(c) 9.50 2.40 0.73 0.17 
\ 
me 
| 
b q_Q-q 
6 Q (14) 
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A significant physical interpretation of the channel-contraction ratio is 
suggested by the last term in Eq. 14. That is, m can be interpreted as 
measure of that part of the total flow which enters from the sides into the 
contracted stream. 

The particular advantage of this concept in connection with the general 
bridge-waterway problem is that the channel-contraction ratio for irregular, 
natural channels can be computed as a ratio of hydraulic conveyances. 
Thus, even for the most extreme combinations of varied depths, alinements, 
and roughnesses in the approach channel, m can be evaluated from 


m=1—£=21- 


Q 
in which K, is the conveyance of that part of the approach channel occupied 
by the discharge g, and Kg is the conveyance of the total section. Con- 
veyance is defined, in terms of Manning’s formula, as 


in which A is the area, and R is the hydraulic radius of the cross sections under 
consideration. 

The General Discharge Equation.—Substituting the three new boundary 
parameters in Eq. 12, the functional relationship for the coefficient of dis- 
charge becomes 
db’ b’ 


c=n(F, m, 


The discharge equation (Eq. 13) is only slightly changed to 


v3, 
Q=Cbys 29 (ah-h, 


To account for irregular cross sections at the constriction, the depth, ys, is de- 
fined as the difference between the water-surface elevation at the downstream 
side of the constriction, adjacent to the contracted stream, and the average 
bottom elevation across the opening. It should be emphasized that y; is 
assumed to be an adequate substitute for y2, the depth of the contracted stream 
at the vena contracta. The difference between these two depths is small for 
all but the larger values of L/b, and even in the latter instance, the difference 
is a function of the variables contained in Eq. 17. Thus, the experimentally 
determined values of the coefficient of discharge always contain the effect of 
this approximation. 


VERIFICATION OF THE APPROXIMATE GENERAL SOLUTION 


Experimental Procedure.—The verification procedure used for the general 
solution involved a systematic investigation of the influence of each of the 
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independent variables contained in Eq. 17. Ideally, this procedure would re- 
quire that the variation of C with each independent variable be measured for 
a full range of each of the remaining variables. An almost unlimited variety 
of conditions would be involved. Actually, considering the minor influence of 
certain variables, the procedure was simplified by analyzing and correlating the 
test results as the experiments were made. The only conditions tested, there- 
fore, were those required to define a series of curves capable of correlating the 
significant variables for a practical range of each. This leads to the procedure 
described previously for the approximate solution of the simple constriction in 
rectangular channels. In this case, however, many more tests were needed to 
develop the average adjustment curves required for the approximate com- 
putation procedure. The range of variables tested and the process of develop- 
ing the empirical relationships are described subsequently. 

Computation of the Friction Head.—It was considered sufficient to investigate 
only three various channel roughnesses in the laboratory. The smoothest 
boundary condition resulted when a cement-asbestos board (smooth surface 
up) was used to create irregular channel cross sections in the flume. Manning’s 
n for this surface, as determined by experiment, was 0.011. Manning’s n for 
the concrete floor of the flume was 0.012, which is not sufficiently different from 
that for the cement-asbestos board to be considered a distinct degree of rough- 
ness. With one layer of 11-gage, 2-in. mesh, chain-link fencing covering the 
floor, the measured n was equal to 0.027. Maximum roughness was created 
by use of two layers of chain-link fencing, for which n was equal to 0.036. The 
net flow area used in computing n for conditions involving the wire fencing was 
the gross area of the section minus the average cross-sectional area occupied 
by the wire. The average depth was computed as the ratio of the net area to 
the top width of the section. 

It is recognized that the principal effect of increased channel roughness is 
to increase the differential head across the constriction. Thus, as roughness 
increases, evaluation of the net velocity head at section 2, 


is increasingly dependent on an accurate evaluation of h;, the loss of head 
caused by friction between sections 1 and 2. For smooth channels and rela- 
tively large values of the channel-contraction ratio, h, is insignificant relative 
to the drop in water surface, Ah. For extremely rough channels, espeically 
those that contain considerable amounts of tree growth and brush, the com- 
putation of friction losses in the accelerated flow approaching the constriction 
is subject to errors of judgment as well as errors in method. 

The procedure used in computing h, in the verification tests was the stand- 
ard method used by the USGS. Briefly, the method consists first of computing 
the average energy gradients at the nominal approach section (1) and at the 
constricted section (2). The energy gradient (slope) is computed as the square 
of the ratio of the discharge to the conveyance of the section under considera- 
tion. The total friction loss is then computed as the head loss in the reach 
between section 1 and the upstream face of the constriction, plus the head loss 
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through the constriction (appreciable only for larger values of L). For the 
upstream reach, the head loss is defined as the geometric mean of the average 
gradients at the end sections times the length of the reach. The loss in the 
constriction is the product of the average gradient at the constricted section 
and the length of the constriction. Acknowledged weaknesses in this method 
of computation apparently did not affect its validity for the conditions tested. 

Location of the Approach Section—The approach section is defined as the 
section which marks the beginning of acceleration of the flow approaching the 
constriction. Since this is not distinctly defined as a line across the flow 
pattern, it is difficult to determine by measurement in the field as well as 
it is in the laboratory. For practical purposes in this case, the beginning of 
acceleration (or beginning of drawdown as it is frequently termed) is of most 
concern at the boundaries. In terms of high-water marks, it is desirable to be 
able to determine the beginning of the drawdown on the upstream face of the 
embankment as well as on the banks of the stream. The error which results 
from inaccurate determination of this location is reflected, in the discharge 
equation, as an error in Ah. When the friction head is small, therefore, this 
error can be avoided by making a generous estimate of the distance from the 
constriction to the approach section. When the friction head is appreciable, 
the assumption of too long an approach section results in a small positive error 
in the computed value of hy. This is seldom of significance in the range of 
conditions considered. 

The location of the approach section is influenced primarily by the Froude 
number and the channel-contraction ratio. Tests were made mainly for the 
purpose of establishing a rough general rule for locating the nominal approach 
section. The tests indicate that a convenient and adequate approximation of 
the length of the approach reach, measured along the center of the channel, is 
a distance equal to the gross width of the opening, b. 

The field conditions for which high-water marks along the upstream face 
of the embankment may ordinarily comprise the best source of data are those 
corresponding to larger channel-contraction ratios. In this case, for values of 
m greater than 0.80, the average water-surface elevation in the approach 
section may be approximated as the elevation of the water at a point a distance 
equal to one half the width (6/2) from the nearest corner of the opening. 

Computation of the Average Kinetic Energy—An additional effect of the 
channel shape and the boundary roughness is reflected in the magnitude of the 
factor a; in the kinetic energy term a x . This term and the term h, in the 
discharge equation, Eq. 17, are generally of approximately equal magnitude 
for the conditions considered. Both quantities are ordinarily small in com- 
parison with the total measured drop, Ah. The coefficient a, was evaluated 
on the basis of detailed velocity surveys for several conditions of flow. It was 
found to vary from 1.06, for relatively large depths of flow in irregular sections, 
to 1.10 for smaller depths in rectangular channels. These average values were 
used in the computations. The error involved in estimating a: for a natural 
ehannel is ordinarily smaller than the error involved in estimating n. 

Evaluation of the Discharge Coefficient—The approximate computation 
procedure described in connection with the verification of the basic analysis will 
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also be used for the expanded analysis. Previously, this procedure necessitated 
the determination of a “standard” empirical relationship between the three 
most critical independent variables contained in the functional equation for C, 
Eq. 14. The variables selected for this purpose are C, m, and L/b. Ignoring 
the ratio y3/b, which was shown by experiment to be insignificant, the standard 
condition is then defined as that for which the remaining ratios are constant 
and equal to the following values: F = 0.5; r/b = 0, e = 1, @ = 0°, and 
W/b = 0. From the experimental data for the standard condition, a family 
of base curves showing the relationship between C, m, and L/b was constructed. 
Individual “adjustment” curves were then developed from systematic in- 
vestigations of the separate influence of the ratios F, r/b, e, ¢, and W/b. In 
each case, as for the basic analysis, the adequacy of the average adjustment 
eurves was then tested throughout a reasonable range of all the related ratios. 

The method of applying the experimentally determined functional re- 
lationships to the computation of the coefficient of discharge is as follows: For 
a given case, C’ is equal to the “standard” value of the coefficient of discharge, 
corresponding to the standard values of F,r/b, e, ¢, and W/b, and the given 
values of m and L/b; kr denotes a coefficient which, multiplied by C’, will 
yield a coefficient of discharge adjusted for the effect of the given Froude 
number, F; k, is a coefficient similar to ke which adjusts the coefficient, of dis- 
charge for the given corner-rounding ratio, r/b; k, represents a coefficient 
similar to ke which adjusts the coefficient of discharge for the given eccen- 
tricity, e; ky is a coefficient similar to kr which adjusts the coefficient of dis- 
charge for the given angularity, ¢; and k,, denotes a coefficient similar to kr 
which adjusts the coefficient of discharge for the given guike-wall ratio, W/b 


(this coefficient cannot exist concurrently with r/b > 0). It follows that the 
adjusted coefficient of discharge can be computed from the equation, 


Cm (rk, by be) C’ (20) 


Effect of the Channel-Contraction Ratio —Experimental verification of the 
eoncept of the channel-contraction ratio in terms of conveyances, Eq. 15, is one 
of the most significant results of this investigation. As applied to channel 
sections of all shapes and roughness patterns, m is a measure of the ratio of 
that. part of the flow which approaches the opening from bank areas upstream 
from the constriction to the total flow which passes through the constriction. 
The coefficient of discharge, C, was previously described as being primarily a 
coefficient of contraction. It is logical, therefore, that with all other inde- 
pendent variables constant, increasing values of m should result in decreasing 
values of C. That is, the degree of ‘‘jet’’ contraction should be proportional 
to the ratio m. Substantiation of this concept establishes m as a general 
boundary-geometry parameter. 

The variation of C with m is shown in Fig. 4(a) for two values of L/b. It 
can be seen that C varies from a maximum value of 1.0 to a minimum of less 
than 0.70. These values, as expected, are in reasonable agreement with values 
of C. for orifices. The range of (L/b)-values for which this relationship was 
tested varied from 0 to 1.0, and the test points shown in Fig. 4(a) are indicative 
of the scatter of the experimental data. 


467-18 


| 

i 

| 
| 


Fig. 4(b), based on a large number of tests similar to those shown in Fig. 
4(a), shows the base curves for the approximate computation procedure. In 
this figure, the coefficient of discharge is indicated by the symbol C’, the 
“standard” value, corresponding to specified constant values of the secondary 
variables, F, r/b, e, ¢, and W/b. The standard values of the secondary 
variables for Fig. 4(b) are the same as in Fig. 4(a). 

It is significant that the definition of the curves shown in Fig. 4(b) was not 
only substantiated but was refined by plotting on the original a great many 
values computed from Eq. 20. Thus, Fig. 4(b) represents the composite of 
experiments made for a wide variety of channel configurations and diferent, 
nonstandard values of the secondary variables. The successful correlation of 
these data constitutes the most critical step in the verification of the approxi- 
mate solution. 

Effect of the Length Ratio —The variation of C with the length ratio L/b is 
indicated by the vertical spread of the curves in Fig. 4(b). For all of the con- 
ditions tested, L was measured as the perpendicular distance from the upstream 
face to the downstream face of the constriction. The width of the opening 
was measured as the minimum clear width between the vertical faces of the 
opening. 

The primary effect of length, as indicated in Fig. 4(b), is to increase the 
coefficient of discharge by decreasing the contraction of the stream as it 
passes through the opening. That is, the contraction is suppressed, even at 
small values of L/b, by the negative pressure gradient which exists between 
the live stream and the eddying fluid which surrounds it. The existence of 
the pressure gradient is-the result of the entrainment of the eddying fluid and 
the reduced opportunity for fluid in the downstream separation zone to replace 
that which is entrained. Thus, as L/b is increased, the area of access for the 
return flow (at the downstream corner of the opening) is increasingly con- 
stricted by the expansion of the stream, and the surface level in zone 3b in 
Fig. 1 tends to drop further below that at the end of the opening, at zone 3a. 
Simultaneously, however, this difference in level is decreased by the energy loss 
which results from the expansion process. As it is the level at zone 3a which is 
used in computing y; in the discharge equation, the effect of the first influence 
is to increase the apparent coefficient of discharge, and the effect of the head loss 
is to decrease the coefficient of discharge. The growth of these conflicting in- 
fluences is indicated by the relative spacing of the curves in Fig. 4(b). 

Effect of the Froude Number—The Froude number can be defined as a 
ratio of a typical unit inertial reaction to a typical unit weight force involved 
in a given fluid motion. Thus, a low value of the Froude number at any 
point in a free-sarface flow pattern ordinarily indicates that gravity, or fluid 
weight, has a large influence on the motion at that point. In this case, the 
prerequisite of tranquil flow limits the Froude number in the contracted stream 
to values less than the critical value. Within this low range, however, larger 
values of F for the contracted stream are usually indicative of a relatively 
large drop in water level, Ah, and a relatively large change in Froude number 
between the approach section and the constricted section for any given bound- 
ary. Furthermore, the effect of gravity on the motion is partly removed by 
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the condition of ‘submergence’ which exists at the downstream side of the 
opening. Thus, the relative influence of the Froude number is proportional 
to Ah, which, for any given boundary geometry, increases with increasing 
values of F for the contracted stream. 

The effect of the Froude number on the coefficient of discharge is indicated 
in Fig. 5, the first of the empirical adjustment curves. Here the value of F is 
the Froude number for the contracted stream and is computed from the 
equation, 


Q 
21 
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Typical data plotted in Fig. 5 demonstrate the apparent lack of correlation 
between kr, the adjustment ratio, and m, the channel-contraction ratio. The 
modifying influence of the remaining independent variables, as determined by 
check tests, is similarly insignificant. The increasing positive slope of the 
adjustment curve with larger values of F is an indication of the increasing in- 
fluence of the Froude number as the flow conditions approach the critical. 

Effect of Entrance Rounding—The most common form of “streamlined” 
entrance for open-channel constrictions is that in which the upstream corners 
are rounded to circular ares. Generally, identical forms are used for both the 
upstream and downstream corners, but the influence of the downstream stream- 
lining is not well established. In this case, the laboratory tests indicated that 
the effect of rounding the downstream corners was negligible. 

The effect of rounding the upstream corners of the constriction was deter- 
mined from tests covering a practical range of the relative radius of curvature, 
r/b. The results of these tests are shown in the form of adjustment curves in 
Fig. 6(a). It was necessary to include a third variable, m, in order to correlate 
the results with the standard values of the remaining independent variables. 

The obvious effect of rounding is to suppress the contraction and increase 
the coefficient of discharge. The curves in Fig. 7 show, in general, that C in- 
creases at a diminishing rate until it reaches a maximum value. This maximum 
value of C corresponds to a value of C, approaching unity. Under these 
limiting conditions, therefore, the coefficient of discharge reflects the boundary 
friction loss. 

Effect of Eccentricity of Opening —The opening is described as being ec- 
centrically located when the length of one of the obstructions which comprise 
the constriction is longer than the other. It would appear that the principal 
effect of eccentricity would be to change the significance of the channel- 
contraction ratio in the functional relationship for C. The experiments show 
that the effect of e on the standard value of C is so small, however, that it can 
be ignored in most cases. The adjustment curve for the effect of eccentricity 
is shown in Fig. 6(d). 

Effect of Angularity——When a constriction is skewed with respect to the 
main thread of the stream, the effective area of the opening isdecreased. It 
follows that the coefficient of discharge is decreased in proportion to the 
magnitude of the skew angle, ¢. In Fig. 7(a), experimental results on the 
effect of @ are plotted as an adjustment curve for the standard value of C. 
The channel-contraction ratio is an important third variable. As indicated 
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in Fig. 7(a), the limiting value of ky, corresponding to m = 0, is equal to the 
cosine of @. For this condition, of course, the channel is not constricted. It 
is apparent from Fig. 7(a) that the effect of angularity is very small for values 
of @ less than 20°. 

Effect of Guide Walls —The influence of diagonal guide walls at the en- 
trance to the constriction is similar to that of streamlining. Only one form of 
guide wall was investigated in the laboratory. This involved plane, vertical 
walls set at an angle of 45° with the face of the constriction. Fig. 7(6) shows 
the- effect of the guide walls on the standard value of the coefficient of dis- 
charge. The shape of the curves in Fig. 7(6) is logically similar to that shown 
by the rounding-ratio curves (Fig. 6(a)). 

Comparison of Adjusted and Observed Values of the Discharge Coefficient.— 
Verification of the approximate solution throughout the complete range of 
variables considered is shown in Fig. 8. In Fig. 8, Cyaj is the coefficient of 
discharge computed from Eq. 20 on the basis of the base coefficient diagram, 
Fig. 4(6), and the adjustment curves, Figs. 5 to 7, inclusive. The “observed” 
value of the coefficient, Co,s, is computed directly from Eq. 17 on the basis of 
the measured values of the primary variables. 

The ratio of C,aj to Cops is a measure of the relative accuracy of the solution. 
This ratio is plotted as a function of C.., in Fig. 8. The indicated accuracy is 
satisfactory within the limits of the conditions tested. These limits are 
tabulated in Fig. 8. 

Illustrative Example-—The recommended procedure for obtaining the ad- 
justed coefficient of discharge from the empirical data is illustrated in the 
following example: These measured data correspond to a test condition ar- 
ranged in a smooth, rectangular laboratory flume: Q = 0.710 cu ft per sec, 
yi = 0.305 ft, ys = 0.255 ft, B = 10.0 ft, b = 2.00 ft, r = 0.083 ft, L = 0.81 
ft, e = 1.00, and ¢ = 0°. Computed data are: Ah = 0.050 ft, hy = 0.001 ft, 

72 
a i = 0.001 ft, m = 0.80, 5 = 0.042, 2 = 0.41, and F = 0.49. For these 
conditions, C.,, computed from Ey. 17 is 0.774. This example involves only 
two adjustment ratios, ke and k,. Taking these from Figs. 5 and 6(a), re- 
spectively, for the conditions indicated, kr = 0.995, and k, = 1.085. The 
corresponding standard value of the discharge coefficient is C’ = 0.705. Thus, 
the adjusted value of the coefficient is Cyaj = ke k, C’ = 0.995 X 1.085 
X 0.705 = 0.761. The relative accuracy of the adjustment procedure in 
this case is indicated by the ratio at = 0.984, which is plotted in Fig. 8 


(“Example”). 
CoNncLUSIONS 
A practicable method of solving the discharge equation for the tranquil 
flow of water through open-channel constrictions is described. The func- 
tional relationship between the coefficient of discharge and the principal 
independent variables is presented in the form of empirical data which were 
the result of a systematic laboratory investigation. Boundary conditions 
considered included vertical constrictive elements, channel cross-sectional 
shapes, and roughness patterns. Within the limits tested, the proposed 
computation procedure yields satisfactory results. 
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Applicability of the proposed solution to bridge waterways in natural 
channels will depend, to large extent, on a critical analysis of the factors 
controlling the flow pattern. For this reason, a significant part of the paper 
is devoted to a detailed discussion of the flow pattern. 
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APPENDIX. NOTATION 


The following letter symbols, adopted for use in the paper and the guidance 
of discussers, conform essentially with American Standard Letter Symbols for 
Hydraulics (ASA-Z10.2-1942), prepared by a committee of the American 
Standards Association with Society representation, and approved by the 
Association in 1942: 


area of channel cross section; 
width of the channel; 
width of the opening; 
coefficient of discharge; 
“standard”’ value of the coefficient of discharge; 
= coefficient of contraction; 
eccentricity ratio; 
= Froude number; 
acceleration caused by gravity; 
piezometric head; 
h. = eddy loss between sections 1 and 2; 
h; = boundary friction loss; 
= conveyance, from Manning’s formula; 
energy loss coefficient ; 
adjustment coefficient relating the adjusted to the standard coeffici- 
ent of discharge; 
k. = adjustment coefficient for the effect of eccentricity ; 
ke = adjustment coefficient for the effect of the Froude number; 
k, = adjustment coefficient for the effect of the corner-rounding 
ratio; 
ke adjustment coefficient for the effect of the guide-wall ratio; 
k, = adjustment coefficient for the effect of angularity ; 
length of the constriction ; 
channel-contraction ratio; 
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\ 
j 
i 
| 
| 
B 
b q 
C 
F 
g 
h 
K 
Ki 
k 
L 
m 


= distance measured normal to the streamlines or the retardation co- 
efficient in the Manning formula; 
= pressure intensity ; 
= total discharge at a channel cross section; 
= unit discharge in a rectangular section, or a part of the total dis- 
charge through a constriction; 
= hydraulic radius; 
= radius of a circular arc; 
= hydraulic slope; 
= distance measured parallel to streamlines ; 
= average velocity; 
v = velocity at a point; 
W = width of the entrance guide walls; 
zx = linear horizontal distance; 
y = depth of flow referred to the average bottom elevation at a cross 
section ; 
elevation referred to an arbitrary horizontal datum; 
= kinetic energy correction factor; 
specific weight of the fluid; 
= density of the fluid; and 
= acute angle between the plane of the constriction and a line normal 
to the thread of the stream. 


| 467-23 . 


Piezometric Head Along 
Outer 


Assumed Separation Boundary 


Plane of Symmetry 


Fic. Constriction 1n A Horizontat CHANNEL 


ELEVATION 


Fic. 2.—Two-Dmmenstonat Constriction wire = 0.5 
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al Stagnation Head, h = ho + ~ | 

Distance Along Center Line or Outer Boundary | L } H | 

! | 

4 

\ 

V =Constant = 

| 

| | 

4s * 

| 

4 

PLAN 

——_—_b— — Average Water Surface 
ah ah 
> 
7 7 
| 


Value of C’ 


Variable Roughness 


(a) 


Fie. 3.— VARIATIONS IN SHAPE OF CHANNELS 


My 


Standard Conditions 


r 


F=0.5 e=l 
$=0° 


45° 


Variable Roughness 


(b) 


(a) C VERSUS m 
x 


Value of 


0 0.50 0.60 0.70 
Value of m 


Fic. 4.—VarraTIon oF Discnuarce Coerrictent C C’ witH m anv L/b 


= (e) 
. 
=. 0.90 
m 
1.00 lose © 
> 
~ 
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Symbol Value of m 
0.20 
x 0.80 


Value of kp 


0.40 0.50 0.60 
Value of the Froude Number 


Fig. 5.—Errect on Discuarce Coerricient oF THE FroupeE NUMBER 


2 
> 


Standard Value 
k,=1.00 
e=1.00 


12 
Value of | (Multiply by 0.01) and e (Multiply by 0.05) 


Fie. 6.—Errect or Corner RounpinG AND ECCENTRICITY 


q 
0.10 020 0.30 0.70 080 090 
{ 
1.20 — 
(a) k, SUS 
0.40 F 
q 
1.05 
GB. 
Standard Value; 0, ky = 1.00 
0.95 | 
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Value of k Py 


8 10 
Value of z (Multiply by 0.01) and m (Multiply by 0.05) 


Fic. 7.—Errect oF ANGULARITY AND Guipr WALLS 


Experimental Limits 
Ratio Range Ratio 


0.20-0.95 e 

0-100 

15-085 F 
-0.13 


Cube 


Fria. 8.—CoMPaARISON OF ADJUSTED AND OBSERVED VALUES OF C 
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uo 


Value of k yy 


1.00 
- 
m=0, ky =Cos 
0.60 
0.30 
| 
1.00 
1.20 — 
om 0 -1.00 
F 0 -0.13 
S t (Manning's) 0.01 1-0.036 
| 
0.60 0.70 a 0.90 1.00 
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PROCEEDINGS-SEPARATES 


The technical papers published in the past year are presented below. Technical-division sponsorship is indicated by an 
abbreviation at the end of each Separate Number, the symbols referring to: Air Transport (AT), City Planning (CP), Con- 
struction (CO), Engineering Mechanics (EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), 
Sanitary Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Water- 
ways (WW) divisions. For titles and order coupons, refer to the appropriate issue of “Civil Engineeriag” or write for a 
cumulative price list. 


VOLUME 79 (1953) 


AUGUST: 230(HY), 231(SA), 282(SA), 233(AT), 234(HW), 235(HW), 237(AT), 238(WW), 239(8A), 240(IR), 241(AT), 242(IR), 
243(ST), 244(ST), 245(ST), 246(ST), 247(SA), 248(SA), 249(ST), 250(EM)®, 251(ST), 252(SA), 253(AT), 254(HY), 255(AT), 
256(ST), 257(SA), 258(EM), 259(WW). 


SEPTEMBER: 260(AT), 261(EM), 262(SM), 263(ST), 264(WW), 265(ST), 266(ST), 267(SA), 268(CO), 269(CO), 270(CO), 271(SU), 
272(SA), 273(PO), 274(HY), 275(WW), 276(HW), 277(SU), 278(SU), 279(SA), 280(IR), 2861(EM), 282(SU), 283(SA), 284(SU), 
285(CP), 286(EM), 287(EM), 288(SA), 289(CO). 

OCTOBER: 290(all Divs), 291(ST)*, 292(EM)®, 293(ST)*, 294(PO)®, 295(HY)*, 296(EM)*, 297(HY)*, 298(ST)*, 299(EM)®, 
300(EM)*, 301(SA)*, 302(SA)*, 303(SA)*, 304(CO)*, 305(SU)*, 396(ST)*, 307(SA)*, 308(PO)*, 309(SA)®, 310(SA)*, 311(SM), 
312(SA)®, 313(ST)*, 314(SA)*, 315(SM)*, 316(AT), 317(AT), 318(WW), 319(IR), 320(HW). 

NOVEMBER: 321(ST), 322(ST), 323(SM), 324(SM), 325(SM), 326(SM), 327(SM), 328(SM), 329(HW), 330(EM)*, 331(EM)*, 
332(EM), 333(EM)°, 334(EM), 335(SA), 336(SA), 337(SA), 338(SA), 339(SA), 340(SA), 341(SA), 342(CO), 343(ST), 344(ST), 
345(ST), 346(IR), 347(IR), 348(CO), 349(ST), 350(HW), 351(HW), 352(SA), 353(SU), 354(HY), 355(PO), 356(CO), 357(HW), 
358(HY). 


DECEMBER: 359(AT), 360(SM), 361(HY), 362(HY), 363(SM), 364(HY), 365(HY), 366(HY), 367(SU)°, 368(WW)°, 369(IR), 
370(AT)©, 371(SM)©, 372(CO)°, 373(ST)°, 374(EM)°, 375(EM), 376(EM), 377(SA)°, 378(PO)°. 


VOLUME 80 (1954) 


JANUARY: 379(SM)°, 380(HY), 381(HY), 382(HY), 383(HY), 384(HY)°, 385(SM), 386(SM), 387(EM), 388(SA), 389(SU)°, 390(HY), 
391(IR)°, 392(SA), 383(SU), 394(AT), 395(SA)°, 396(EM)°, 397(ST)°. 


FEBRUARY: 398(IR)4, 400(CO)4, 401(sM)°, 402(AT)4, 403(AT)4, 404(IR) 4, 405(PO)4, 406(AT)4, 407(sU)4, 408(SU)4, 
400(ww)4, 410(A7T)4, 411(Sa)4, 412(P0)4, 


MARCH: 414(ww)4, 415(sU)4, 416(sm)4, 417(sM)4, 418(AT)4, 419(8a)4, 420(SA)4, 421(AT)4, 422(SA)4, 423(CP)4, 424(AT)4, 
425(sM)4, 426(IR)4, 427(ww)4. 


APRIL: 428(HY)°, 429(EM)°, 430(ST), 431(HY), 432(HY), 433(HY), 434(ST). 
MAY: 435(SM), 436(CP)°, 437(HY)°, 438(HY), 439(HY), 440(ST), 441(ST), 442(SA), 443(SA). 


JUNE: 444(SM)®, 445(SM)®, 446(ST)®, 447(ST)®, 448(ST)®, 449(ST)®, 450(ST)®, 451(ST)®, 452(SA)®, 453(SA)®, 454(SA)®, 
455(SA)®, 456(SM)®. 


JULY: 457(AT), 458(AT), 459(AT)°, 460(IR), 461 (IR), 462(IR), 463(IR)°, 464(PO), 465(PO)°. 


AUGUST: 466(HY), 467(HY), 468(ST), 469(ST), 470(ST), 471(SA), 472(SA), 473(SA), 474(SA), 475(SM), 476(SM), 477(SM), 
478(SM)°, 479(HY)°, 480(ST)©, 481(SA)°, 482(HY), 483(HY). 


a. Presented at the New York (N.Y.) Convention of the Society in October, 1953. 

b. Beginning with “Proceedings-Separate No. 290,” published in October, 1953, an automatic distribution of papers was in- 
augurated, as outlined in “Civil Engineering,” June, 1953, page 66. 

c. Discussion of several papers, grouped by Divisions. 

d, Presented at the Atlanta (Ga.) Convention of the Society in February, 1954. 

e. Presented at the Atlantic City (N.J.) Convention in June, 1954. 
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